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Abstract. In this paper we investigate the idea of Hanasz 
& Lesch 1993 that the galactic dynamo effect is due to the 
Parker instability of magnetic flux tubes. In addition to 
the former approach, we take into account more general 
physical conditions in this paper, by incorporating cos- 
mic rays and differential forces due to the axisymmetric 
differential rotation and the density waves as well. 

We present the theory of slender magnetic flux tube 
dynamics in the thin flux tube approximation and the 
Lagrange description. This is the application of the for- 
malism obtained for solar magnetic flux tubes by Spruit 
(1981), to the galactic conditions. We perform a linear sta- 
bility analysis for the Parker-shearing instability of mag- 
netic flux tubes in galactic discs and then calculate the 
dynamo coefficients. 

We present a number of new effects which are very 
essential for cosmological and contemporary evolution of 
galactic magnetic fields. First of all we demonstrate that 
a very strong dynamo a-effect is possible in the limit of 
weak magnetic fields in presence of cosmic rays. Second, 
we show that the differential force resulting from axisym- 
metric differential rotation and the linear density waves 
causes that the a-effect is essentially magnified in galactic 
arms and switched off in the intcrarm regions. Moreover, 
we predict a non-uniform magnetic field in spiral arms and 
well aligned one in interarm regions. These properties are 
well confirmed by recent observational results by Beck & 
Hoernes (1996) 
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1. Introduction 

The origin of galactic magnetic fields is a long standing 
problem in theoretical astrophysics. Virtually all spiral 
galaxies have magnetic fields of a few /iG which are spa- 
tially coherent over several parsecs (see Kronberg 1994; 
Beck et al. 1996 for recent reviews). This implies that the 
average magnetic energy density is comparable to the av- 
erage energy densities of cosmic rays and the interstellar 
medium, which in turn suggests that galactic magnetic 
fields reached a state of saturation some unknown time in 
the past. In face of the Faraday rotation measurements of 
Wolfe et al. (1992), which indicate /iG fields in damped 
Lymana clouds at redshifts of 2 (the universe was about 
1-2 Gigayears old), the saturation has to be very rapid. 
Since galaxies are so large (and damped Lymana clouds 
are supposed to be progenitors of disc galaxies), and their 
constituent gases are highly conducting, the dissipative 
time scale for these fields is considerably longer than a 
Hubble time. Either the magnetic fields are simply the re- 
sult of some pre-existing cosmological field (Kulsrud 1990) 
or they are spectacular examples of fast dynamos, in which 
rapid magnetic reconnection takes over the role of the tra- 
ditional eddy diffusivity (Parker 1992) 

Nowadays the regular magnetic fields in the discs of 
spiral galaxies are usually considered to be the result of 
large-scale dynamo action, involving a collective induc- 
tion effect of turbulence (a-effect) and differential rota- 
tion. Both mechanisms have been described in detail in a 
huge number of articles (e.g. Kronberg 1994; Beck et al. 
1996 and references therein). A small seed field is exponen- 
tially amplified by the interplay of turbulent motions per- 
pendicular to the disc generating a radial magnetic field 
component By from an azimuthal component B^. Differ- 
ential rotation in the galactic disc then closes the loop by 
generating from Br. The magnetic field grows exponen- 
tially on time scales of the order of a Gigayear (Camenzind 
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and Lesch 1994). The growth time raises the problem of 

magnetic field amplification at high redshifts. The typical 
seed field strengths is of the order of 10~^* G and the typ- 
ical present magnetic field strength is of the order of 10~^ 
G. The e- folding times of the dynamo magnification pro- 
cess have to be of the order of a few hundreds of Myr or 
even a few tens of Myr for the total time of growth 10 Gyr 
or 1 Gyr respectively. Such a short c-folding time seems to 
be unavailable within the classical dynamo models (Lesch 
and Chiba 1995). 

Whereas the differential rotation of galactic disc is a 
clearly observed feature of spiral galaxies, there is no gen- 
eral agreement on the source or value of a. The helicity a 
is connected with the interstellar turbulence, which may 
be fed by several sources (warm gas, hot bubbles, cold fil- 
aments, supernova explosions, cloud-cloud collisions, OB- 
associations, the galactic fountain, etc.). These violent 
motions agitate the interstellar gases and drive turbulence. 
The classical a-effect measures the field-aligned electro- 
motive force resulting from magnetic field lines twisted 
by the turbulence. Originally the dynamics of these field 
lines is governed by external turbulent motions like super- 
nova explosions (Ruzmaikin, Sokoloff and Shukurov 1988; 
Ferriere 1993). However, attempts to compute the sev- 
eral components of the helicity tensor ay resulting from 
exploding supernovas or supcrbubbles resulted in values 
considerably to small to allow for fast and efficient dy- 
namo action (Kaisig et al. 1993). Nevertheless, the most 
recent calculation of the superbubble contribution to the 
helicity tensor by Ferriere (1996) makes the superbubble 
model more attractive. Beside that time scale problem, 
the recent observations by Beck and Hoernes (1996) raise 
the question what kind of global process in a galactic disc 
can explain the "magnetic arms" in NGG 6946. They re- 
port the detection of two highly polarized magnetic fea- 
tures which are 0.5 - 1 kpc wide and about 12 kpc long 
and have greater symmetry than the optical arms. These 
features cannot be explained by standard dynamo models. 

The usual dynamo models treat the magnetic field as a 
superposition of a homogeneous mean field and a fluctuat- 
ing, turbulent field component, which is responsible for the 
turbulent diffusion and the helicity, as well. We considered 
in a first paper (Hanasz and Lesch 1993, hereafter HL'93) 
the Parker instability of magnetic flux tubes (Parker 1955, 
1966, 1967a,b) as the source for the helicity, taking into ac- 
count the observed complicated polarized filaments, which 
resemble the complex structure of the interstellar gas ve- 
locities and densities. The concept of flux-tube dynamo is 
well known for stars (Schiissler 1993), but has never been 
developed for galaxies. 

The Parker instability occurs in a composition of gas, 
magnetic field and cosmic rays placed in a gravitational 
field, if the magnetic field and the cosmic rays contribute 
to a pressure equilibrium of the system. The given pressure 
equilibrium of the flux tube implies that the density of the 
composition of gas, magnetic field and cosmic radiation is 



smaller than the density of the surrounding gas, which in 
the presence of a gravitational field leads to buoyant field 
lines. This situation is clearly realized in galactic halos. If 
the Parker instability operates, the gas together with the 
frozen-in magnetic field and cosmic radiation captured by 
the magnetic field start to ascend, forming bumps on the 
initially straight horizontal field lines. The portion of gas 
in the ascending area starts to sink along the field lines, 
forming the condensations of the gas in the valleys and 
rarefaction in summits. 

Taking into account the dragging effect of turbulence 
in the surrounding medium we could show that the Parker 
instability provides a helicity of about a ~ 0.6 km s~^ 
and that the turbulent diffusion coefficient is reduced by 
the effect of aerodynamic drag force, which causes an in- 
crease in dynamo efficiency. This means that a dynamo 
driven by Parker instability saturates faster than a clas- 
sical dynamo. It also means that the magnetic field in a 
galaxy will always stay in the saturated state, regulated by 
magnetic flux loss via the Parker instability. Similar ideas 
for a fast dynamo were proposed by Parker (1992), who 
suggests a mechanism with rising and inflating magnetic 
lobes driven by cosmic rays. 

In our first article we did not consider the effects of the 
cosmic rays and the underlying axisymmetric differentially 
rotating disc, including non- axisymmetric density waves. 
It is the aim of this contribution to show the results of an 
extended version of a flux tube dynamo, including cosmic 
rays in a differentially rotating disc with density waves. 

The influence of the cosmic rays is obvious in face of 
galaxy evolution. Simply the fact that typically 90% of 
the gas in a galaxy has been transformed into stars means 
that star formation was very efficient in the early phases 
of evolution. For the magnetic field evolution in the flux 
tube model this means that the Parker instability was ef- 
fectively driven by the large cosmic ray flux, since a strong 
star formation is always accompanied by many supernova 
explosions, i.e. by efficient particle acceleration. Since the 
cosmic ray pressure can be assumed to be dominant in 
star forming galaxies (see also M82 (Reuter et al. 1992)), 
it regulates the efficiency of the Parker instability. 

The role of differentially rotating disc on the Parker 
instability is rather complicated. The differential forces 
introduce a new level of complexity to the problem. Ac- 
cording to Foglizzo & Tagger (1994, 1995, hereafter FT'94 
and FT'95) the action of differential rotation depends on 
the magnitude of shear but also on the radial and vertical 
wavenumbers of the perturbation. The differential rota- 
tion can transiently stabilize mode of long vertical wave- 
length, while the waves with short vertical wavelength 
are subject to both the Parker and the transient shear- 
ing instabilities. In the case of general 3D distribution of 
magnetic field the linear coupling of different modes takes 
place depending on the actual range of parameters in- 
cluding the components of the wavenumber. These results 
are not easily portable to the flux tube case because of 
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the geometrical differences. Nevertheless, we confirm the 

close "cooperation" of the buoyancy and shearing forces 
and following FT'95 apply the unified name: "the Parker- 
shearing instability". It is however more appropriate to 
use the name "Parker instability" in some contexts. 

The magnetic shearing (Balbus-Hawley) instability 
has been discussed in connection with the buoyancy and 
dynamo action by Brandenburg et al. (1995). The authors 
simulate the non linear evolution of three-dimensional 
magnetized Keplerian shear flow. Their system acts like 
dynamo that generates its own turbulence via the men- 
tioned instabilities. They observe signatures of buoyancy 
as well as the tendency of magnetic field to form intermit- 
tent structures. 

The plan of the paper is as follows: In Section 2 we 
present a general formalism for the description of the flux 
tube dynamics in galactic discs using the thin flux tube 
approximation. In Section 3 we perform a linear stability 
analysis and discuss the relations between the flux tube 
and the continuous approaches to the Parker-shearing in- 
stability. We calculate first the dynamo coefficients of the 
galactic fiux tube dynamo for the case without shear. Then 
we take into account differential forces due to the galac- 
tic axisymmetric shear and the non axisymmetric density 
waves as well. In Section 4 we demonstrate relations be- 
tween our theoretical and the recent observational results 
and propose a new dynamo model based on the Parker- 
shearing instability in galactic discs with density waves. 



of time t and a "spacelike" parameter m - the mass mea- 
sured along the tube. The lagrangian vectors 



dx 
'dt 

dx 
dm 



dy\ f dz 
dm / ' \ dm 



(2) 
(3) 



are respectively a velocity vector and a tangent vector. 
Let us specify the parameter m by the requirement that 
the mass flow across points m = const vanishes, what 
means that we are going to use the lagrangian description. 
Having a freedom to rescale the parameter m we can chose 
the normalization 



dm = ^ds, 



(4) 



where ^ and ds are respectively a mass per unit length 
and a length of an element of the flux tube. Since 



(5) 




then 
\L\= 



(6) 



2. Magnetic flux tubes in galactic discs — basic 
equations and properties 

We are going to apply the thin, slender flux tube approxi- 
mation which according to Spruit (1981), Spruit and Van 
Ballegooijen (1982) means that the tube is assumed to 
be thin compared with both the scale height H of verti- 
cal density stratification and with the scale L along the 
tube on which the direction of its path changes. In the 
foregoing considerations we shall follow the formulation 
of Moreno-Insertis (1986) with some modifications related 
to the specifics of the galactic disc. We will neglect the 
details of the shape of the cross section of the flux tube 
and assume that it is circular. We will also assume that 
the tube is not twisted. The assumptions of the concept 
of the the thin, slender flux tube are extensively discussed 
by Schramkowski and Achterberg (1993). 

2.1. Geometry of the thin flux tube 

Let us describe the shape of flux tube in the 3-D volume 
and the cartesian reference frame rotating with the disc, 
by means of 3 functions 



The unit tangent vector I is given by 
l = ^L. 

Now let us define the curvature vector 



K 



m 

ds 



Then the curvature radius is given by 



1 
R 



= K=\K\=^ 



di 

dm 



and 



n 



K 

'k 



(7) 



(8) 



(9) 



(10) 



is a unit normal vector (principal normal) parallel to the 
curvature vector K. In the following considerations we 
will need to decompose vectors as eg. gravitational accel- 
eration g into components parallel and perpendicular to 
the unit tangent vector I 



X = x{t,m), 
y = y{t,m), 
z = z{t,m) 



(1) 9\\ = {9-1)1, 

9± = (1x9) X I- 



(11) 
(12) 
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2.2. MHD equations for a thin flux tube 



If we assume for simplicity that 5 is a constant, then 



In order to describe the dynamics of magnetic flux tubes 

in a rotating reference frame of galactic discs we shall start 
with MHD equations in the following form 



dp 
dt 



+ 6iy{pv) = 0, 



(13) 



Peiz) = /9e(0)exp 



\z\ 

ll 



(20) 



where H is the scale height of the disc given by Parker 
(1966) 



p[% + {v V). 

-Vp + pg+ -^3 X B -2pnxv + Fdis, 
47r 



dB { 1 

— — = rotfv X B) - rot ( - — rotS 

divB = 0, 



(14) 



(15) 
(16) 



where we have used traditional assignments for physical 
variables. The terms on the rhs. of the equation of motion 
(14) related to rotation need an additional explanation. 
The —2/3$! X t> is the Coriolis force. The additional term 
-F'diff is the differential force due to counteracting effects 
of centrifugal force and the radial component of gravita- 
tion in the shearing sheet approximation (see eg. FT'94 
and references therein). The differential force will be in- 
troduced in the linearized form in the section 3. The role 
of the differential force due to the galactic dynamics is 
discussed in more details by FT'94. In addition to the 
above system of MHD equations we should specify ad- 
ditional relations which are typical for the galactic disc. 
Our intention is to define a general setup for the Parker- 
shearing instability of magnetic flux tubes which is as close 
as possible to the traditional continuous 3D distribution 
of galactic magnetic fields discussed recently in FT'94 and 
FT'95. This will allow us to show that the flux tube sta- 
bility properties are physically very similar to stability 
properties of the mentioned continuous distribution of in- 
terstellar magnetic fields. 

Following Parker (1966, 1967a, b) we shall assume the 
following isothermal equation of state for interstellar gas 



Vg 



u'^p, 



(17) 



with the sound speed of gas u independent of height in 
the galactic disc. 

The magnetic field and cosmic ray pressures in the un- 
perturbed state are traditionally assumed to be propor- 
tional to the gas pressure 



Pmag — ^Pg^ Per — PPg 



(18) 



where a and (3 are constants order of 1. 

The vertical equilibrium of the composition of gas. 
magnetic field and cosmic rays in the vertical gravitational 
field g = —sign{z)\g\ez is described by the equation 



^ {Pg + Pmag + Per) = p{z)gz 



(19) 



H ={l + a + P) — . 



(21) 



In the following we shall define our model of galactic 
flux tubes, which is in fact slightly different with respect 
to that proposed in HL'93. The present model assumes 
that the flux tubes 

1. are composed of magnetic field, ionized gas and cosmic 

rays, 

2. can execute motions which are almost independent on 
the ambient flux tubes (eg. a particular flux tube can 
rise due to the Parker instability, while the neighbor- 
ing flux tubes are temporarily in a kind of dynamical 
equilibrium) . 

3. are distributed in space with the filling factor close to 
1, 

4. are initially azimuthal. This simplifying assumption 
will allow us to construct the equilibrium states neces- 
sary for the linear stability analysis. 

5. Collisions of the flux tubes are temporarily ignored. 
They can lead to the magnetic reconnection processes 
and some reorganization of the magnetic field struc- 
ture, but this topic will be discussed elsewhere. 

Then, we can introduce a division between quantities spe- 
cific to the flux tube interior (z) and its exterior (e). Taking 
into account various component of pressure we can postu- 
late the following magneto-hydrostatic balance condition 



Pig ~t~ Pim ~t~ Pier — Peg ~t~ Pem ~t~ Peer i 



(22) 



stating that the total internal and external pressures are 
the same at the cylindrical flux tube boundary and the 
aerodynamic drag force. The aerodynamic drag force will 
be considered for completeness in our set of equations, 
however it will not be taken into account in the solutions 
of the present paper. We shall operate within a frame of 
linear stability analysis, so the aerodynamic drag force, 
quadratic in velocities, will be neglected. For the drag force 
Fd wc will adopt the formula 



' D 



CoPe 



(23) 



where F'£) is the drag force per unit length (Schiissler 
1977, SteUa & Rosner 1984). Cd is a drag coefficient 
which, over a wide range of physical regimes character- 
ized by subsonic motions and for Reynolds numbers lying 
between 30 and 10^ (Goldstein 1938) is estimated to be 
~ 1 - 10. 
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Now we are going to reduce the equations to the conve- 
nient form for describing the flux tube dynamics. Combin- 
ing the continuity equation (13), the induction equation 
(15) and assuming infinite electrical conductivity one ob- 
tains 



D fB 

Dt 1^7 



B 



• V U 



(24) 



Substituting i? = $/S and p = ^/S, where $ is the flux 
across the flux tube, S is the cross section area, we can 
write 



which leads finally to the equation 
DL , dv 



(25) 



(26) 



In the absence of CorioUs force the equation of motion 
(14) can be decomposed into two components: parallel and 
perpendicular to the tube axis 



Dt J II Pi dm " Pi 



(27) 



Dv 
'Dt 



^K+P^g^ + ^ + ^^, (28) 
47rpi Pi Pi Pi 



where F^iff is the differential force due to the galactic 
dynamics, which will be introduced in more details later 
on. The presence of Coriolis acceleration 



-2n X V 



couples the parallel and perpendicular components of the 
equation of motion. 

Let us assume that the flux tubes are aligned horizon- 
tally in the initial state and assign Pcr{z0) = f3u'^pe{zo) 
to the cosmic ray pressure specific for the height zo in 
the galactic disc. Following Shu (1974) we assume that 
the pressure gradient of the cosmic ray gas is orthogonal 
to the magnetic field lines. Then, during evolution of the 
particular tube the internal cosmic ray pressure remains 
constant and equal to Pcr(zO). Thus, we can rewrite the 
magneto-hydrostatic balance condition (22) as follows 



u^Pi{t,m) + 



Stt 



+ pu pe{zo) 



(30) 



u^Pe{z{t,m)){l + a + l3). 



Due to the magnetic flux conservation we can rewrite the 
hydrostatic balance condition in the form 



1 



Pi = Pe 



u SttE^ 



(31) 



where we denoted 

Pe = Pe{z){\ -I- a + /?) - l3pe{zo) 



(32) 



and = i?fo^o "^it^ Bfg = Snau'^ pe{zo) . Since ^ = piT, 
then the equation relating the two unknowns ^ and S is 



PeS^ - - ape{zo)J:l = 0. 
We find that 

_ ^ + ^e + '^ape{zo)Pe^ 

~ We ■ 



(33) 



(34) 



The final form of the system of equations is as follows 
Dx 



Dt ' 
DL _ dv 

'Dt ~ dm' 
Dv 



dpi 



— = -u^E^l+g 



Dt 



dm 



(35) 
(36) 

^K+Pl^g^ (37) 



.2nxv + ^ + ^. 

Pi Pi 



2.3. Basic parameters and physical units 

Let us now specify the values of basic physical parameters 
involved in our considerations. Following Parker (1979) p. 
806 we shall use parameters typical for the Milky Way as 

an example and take 



3 2:: 2 • 10"^cm s"^. 



(29) p ~ 1.6 • 10"2^g cm 



(38) 
(39) 



which is equivalent to 1 H-atom cm ^. The conventional 
value of the magnetic field strength is 



B ~ 3 • 10-^G 



(40) 



which leads to a magnetic pressure Pm = B'^ /Sir = 0.4 • 
lO^^^dyn cm~^. For comparison the typical value of the 

cosmic ray pressure is 



Per ^ 0.5 • 10"^^dyn cm"^ 



and a gas pressure 

Pg = pu^ ~ 0.4 • 10"^^dyn cm'^ 



(41) 



(42) 



where u ~ 7km s~^ is the thermal velocity of gas. It is 
worthwhile to notice that Parker considers the gas pres- 
sure due to the both thermal and turbulent gas motions. 
Since the above estimations of gas, magnetic field and cos- 
mic ray pressures give almost identical contributions of all 
the 3 components then 



a ~ /? ~ 1 



(43) 
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and the resulting vertical scale height is 

iJ ~ 2.3 • lO^pc. (44) 

In the following considerations we shall use convenient 
units: 1 pc as a unit of distance and 1 Myr as a unit of -^j^^iiici^ implies that 
time. The unit of velocity is 

1 pc Myr"^ ~ 0.98 km s'^ 

and the unit of acceleration is 

1 pc Myr-2 ~ 3.09 • lO"^"* km s~^ 

In these units 

g ~ 0.65 pc Myr"^ 

and 

17-3.2- 10-2 Myr-^ 



(45) 
(46) 
(47) 
(48) 



Using the equation (36), one can find that: 

d dmi 
dt dmo 



mi = 0. 



(60) 



(61) 



The mass per unit length of the tube in the linear approx- 
imation is 



e = eo 1-6 



dxi 
^ dmo 



Now we calculate pn in the following sequence 



dxi 
dmo ' 



3. The linear stability analysis 

3.1. Equations 

With the aim of performing the linear stability analysis 
we decompose all the quantities in the 0-th and 1-st order 
parts. Let us start with 



m = mo + TOi (too , t) 
which leads to 



d 

dm 



= 1- 



dmi{mo,t)\ d 



dmo 



dmo 



(49) 



(50) 



For the flux tube which is extended along the .T-coordinatc 
in the unperturbed state, the coordinates of a given flux 
tube element in the linear approximation are 



mo , 



(51) 
(52) 



y = yi, 

z = zi. (53) 

The components of the tangent vector in the linear ap- 
proximation are 

dxi 



LxO — 
L 
L 



-'yo 



1 

1^ 

dyi 
dmo 

dzi 

dmo 



1 - 



dmi\ 
dmo 



) dmo ' 



(54) 
(55) 
(56) 



The modulus of the tangent vector, the unit tangent vector 
and the curvature vector are respectively 

dmi 



L\ = 
I = 
K = 



dmo 



° dmo ) 



dyi dzi 
l>?o-5 — ,?o 



dmo dmo 
d'^yi d'^zi 



dmo^ ' dmo^ 



(57) 
(58) 
(59) 



Pel = -PeUoj-o^l) 



Si 



Pil 



-'O 



2a + 1 
2a 



Ci , 1.91 



^iiki 

2a -M So 2a -Fl So 



zi 



(62) 

(63) 
(64) 
(65) 
(66) 



Let us derive some coefficients appearing in the equation 
of motion (37) 



Pt - Pc 



Pi 

\4npi 



K 



-1 |g| 
2a + 1 u2 



^1 zi 



2a 6 
2a + Uo ' 



(67) 



(68) 



where v\ = 2av? . The equation of tube motion (37) after 
linearization is 



dvi 
'df 



d'^xi (1 + a + /3) dzi 



+ 



l + 2a\ ds^ 
1 + a + (3 f v\ dxi 



H 



ds 



lo 



(69) 



l-h2a 



Hds -(-^^H^^ 



+ v\Ki -2nxvi+ ttdiff, 

where s is the length parameter of the tube defined by (4) 
and 



—2fl X vi = 



9 = [0,0,9], 

10 = [1,0,0], 
d% dhi 



' ds 



2 

dyi 



ds' 



(70) 
(71) 

(72) 
(73) 



The linear eflfect of an axisymmetric differential rota- 
tion can be incorporated in a way similar to FT'94. The 
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differential force is radial in direction and is proportional 
to the radial displacement in the shearing sheet approx- 
imation. There is an obvious difference between the 3D 
continuous case and the flux tube approach. In the con- 
tinuous case the differential rotation influences also the 
radial wavenumber which varies linearly as a function of 
time (FT'94). In the flux tube approach the wavevector of 
the perturbation is always parallel to the tangent vector 
and this effect is absent since we consider (until now) only 
azimuthal flux tubes which are initially localized radially 
and vertically. 

The acceleration operating on an element of the flux 
tube in the shearing sheet approximation is 



Cldiff = 

where 



Ro dn 
~2dR 



(74) 



(75) 



is the Oort constant {A/Q = —1/2 for flat rotation and 
—3/4 for the Keplerian rotation). 

3.2. The linear .solution 

In a case of vanishing rotation and shear we assume the 
harmonic perturbation on a magnetic flux tube placed in 
the coordinate system (x, z) in the form 



Xi 













exp i{ks — ujt) + c.c. 



(76) 



where ,s is a length parameter measuring position along the 
flux tube, k is the wavenumber, x and z are the azimuthal 
and vertical coordinates respectively. The linearized equa- 
tions describing evolution of harmonic perturbations on 
the flux tube lead to the linear dispersion relation 



det 



l+2a 
l+a+f3 v'a 
l+2a 



l+2a 

l+2a 



H 



same as for no rotation. In the continuous case of FT'94, 

the marginal instability point is the same as for no ro- 
tation {kp) for the radial wavenumber /c^ ^ oo, and is 
smaller [kg < kp) for kr = 0. But in the limit kr — > oo, 
the dispersion relation of FT'94 is independent of the rota- 
tion (A(a;'^) = 0), and of second order in u"^ (their eq. Al). 
The main difference with our approach is that the optimal 
wavenumber k leading to the maximum growth rate is in- 
dependent of the rotation in the continuous case, whereas 
it depends on the rotation in the flux tube approach ac- 
cording to our Fig. 1. This is mathematically different, 
but physically acceptable, since the trends observed in our 
Fig. 1 and in Fig. 1. of FT'94 are very similar. Thus, we 
can conclude that the slender flux tube model, in addi- 
tion to its observational justification, has the advantage 
to reduce the fully three-dimensional MHD equations to a 
curvilinear one-dimensional formalism. Similar conclusion 
has been already made by Schramkowski and Torkelson 
(1996) in the context of accretion discs. 

In the presence of nonvanishing galactic rotation one 
should take into account also the radial coordinate y. In 
this case the linear harmonic perturbations are assumed 
to be 



= 0,(77) 



where ?i and va are the thermal and Alfven speed re- 
spectively, and H = {1 + a + (5)u'^ / \ g \ is the vertical 
scaleheight. After the substitution of the dimensionless: 
azimuthal wavenumber k = kH and complex frequency 
a) = u)H/va, the above dispersion relation reduces exactly 
to the form (Al) in FT'94, which represents a case of con- 
tinuous 3D distribution of magnetic field with vanishing 
vertical wavenumber kz = and the radial wavenumber 
kr = oo. This means that our flux tube approach is analo- 
gous to the above limiting case of continuous distribution 
of magnetic field. 

Nevertheless, as soon as we take into account rotation 
(see below), the analogy becomes less clear. Our dispersion 
relation becomes of third order in oj"^, and the marginal 
stability point kmarg (which can be determined after the 
substitution of w = to the dispersion relation) is the 



Xl 






yi 




Yi 


.-^1 . 







exp i{ks — uit) + c.c. 



We can write the dispersion relation in the form 



(78) 



det 



\-\-2a 

l+a+P v" 



2iQ.u) 



l+a+/3 v\ -J 
l+2a H 



1 + 

= 



2q 



"-^i-tk) 



iAn -uj^ 

(a-/3)(l+a+/3) «^ 
l+2a 







(79) 



If w is a solution of the above equation representing the 
Parker mode for given k, then eigenvectors are given by 
the following relations between constants Xi, Yi and Zi 



Xi_ 
Zi 

Yi_ 

Zi 



(a-£)(l+a+^_ui 
l+2a ff2 



l+2a H y '^ l 



2iU,L0 



Xl 



vik^ + AAn - u;^ Zi 



(80) 
(81) 



In the Lagrangian description, the dynamo coefficients 
ad and % can be calculated using the method applied by 
Ferriz-Mass et. al (1994). The expression for aa is 



(82) 



where b is the first order perturbation of magnetic field, 
V = r is the velocity of a lagrangian element of the fiux 
tube and r = dr/ds is the tangent vector. The angle 
braces (• • •) stand for space averaging, which in our case 
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is related to integration over the parameter s. Similarly, 
the expression for diffusivity is 

Vd = {r-r). (83) 
Assuming only vertical diffusion for simplicity we get 

Vd = ihr,). (84) 

This simplification is valid if horizontal displacements of 
the flux tube are smaller than the vertical ones. We shall 
indicate a case which docs not fulfill this requirement. 
After substitution of the linear solution one obtains 



ad = UikZi{tf Im. 

Vd = Y^iitf 



(85) 

(86) 



where Zi(t) = ZiCxp^LUit) is the amplitude of vertical 
displacement of the flux tube. The neglected horizontal 
diffusivity = ^Yi{t)^ can be easily estimated basing 
on vertical diffusivity since rjdh = ildv\Yi/Zi\'^ . 

The magnetic Reynolds numbers and the dynamo 
number are defined as 



R„ = 



Vd 

R - 

Hui = 

m 

D = RaRoj 



(87) 

(88) 
(89) 



The coefficients given by the relations (85)-(89) de- 
pend on all the involved parameters, especially on the 
vertical density gradient and the magnitude of shear im- 
plicitly, via the solutions of the dispersion relation (79), 
oj = uj{k) and the relations (80) and (81) between com- 
ponents of the eigenvectors. The most relevant quantities 
contributing to ad and f]d will be presented in a series of 
forthcoming figures. 

S.3. The numerical results 

Fig. 1 shows the growth rate vs. wavenumber for a = /3 = 
1, and other parameters are described in the subsection 
2.3. The full line is for the case without rotation, the 
dotted line is for the rotation = ^Iq = 10~^^s~^ = 
0.03Myr~^, the dashed line is for O = 2Qg, and dashed- 
dotted line is for fl = ^flc- As is well known, increasing 
the rotation frequency f2 dimishes the growth rate. 

It is worth noting, however that the ratio of amplitudes 
Yi/Zi of radial to vertical displacements which measures 
the magnitude of cyclonic deformation of the flux tube 
increases with rotation fl (see Fig. 2, where dotted, dashed 
and dotted-dashed curves represent the same values of Q, 
as in Fig. 1). This ratio contributes to the ad coefficient 
in (85). In addition, in the limit of A; — > 0, Yi/Zi — > 1 for 
all values of the rotation frequency. 
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Fig. 1. The dependence of growth rate uji on wavenumber k 
for different values of the galactic rotation frequency f2 
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Fig. 2. The dependence of the ratio of amplitudes Yi/Zi on k 
for the same values of values of Q as in Fig. 1. 



The dependence of ad on k is shown in Fig. 3 for 

= Via/ 2 (continuous line), Vl = flo (dotted line), 
= 2Qg (dashed line) and f2 = 3Qg (dashed-dotted hue). 
The ad coefficient is computed with the assumption that 
the uppermost part of the tube has passed the vertical 
distance equal to the vertical scale height H: Z(t) = H 
for each value of k (for a = p = l, u = 7 km s~^ 
and \ g \= 2 ■ 10~^ cm = 0.65 pc Myr^^ we have 
H = 226 pc). This calculation of ad following Ferriz-Mas 
et al. (1994) is intended only to give a first insight about 
the dependence of ad on k and f7. It hides, however our 
ignorance about: (1) the initial amplitude of perturbations 
and (2) the final state reached by a single flux tube. Thus, 
within the linear theory we are not able to average over a 
statistical ensemble of flux tubes. Fig. 3 shows that the ob- 
served value of the galactic rotation in the solar neighbor- 
hood approximately maximizes the magnitude of the ad 
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coefRcient. This property could be quite important since 

many qualitative approaches adopt ad proportional to fi 
(see eg. Ruzmaikin et al. 1988). What we observe in our 
model is in contradiction with this rule. 



3r 




0.000 0.002 0.004 0.006 0.008 0.010 
k [pc '] 

Fig. 3. The dependence of the dynamo-Qd coefRcient on k for 
a fixed vertical displacement Z\ = H and fi = 1/2, 1, 2, 3 • fie- 

It is worthwhile to notice that the sign of is positive 
for positive z coordinate, as it follows from the nature of 
Coriolis force which tends to slow down the rotation of 
rising and expanding volumes. The flow of gas directed 
outward of the top of the rising flux tube can be recog- 
nized as expanding in the direction along the magnetic 
field lines. 

This feature of our model is in contrast with the results 
of Brandenburg et al. (1995). Strictly speaking they dis- 
cuss two quantities: the mean hclicity < v-votv > and the 
ocdyn defined as a coefRcient relating the azimuthal com- 
ponent of electromotive force with the azimuthal mean 
magnetic field. While the first quantity has a proper sign 
('-' in their convention of assignments, which is consistent 
with the positive sign of our a^), the negative sign a^yn 
coefficient seems to be in conflict with the positive sign of 
helicity. 

We explain that the coefficient in our model is an 
equivalent of the hclicity < v ■ rotw > used by Brandenburg 
et al. (1995). On the other hand we are dealing with a 
simple model in which the dynamo properties are derived 
from the dynamics of a single flux tube. Studies of the 
dynamics of an ensemble of flux tubes would help to find 
a reason of the discrepancy. Moreover, we remain in the 
linear regime in the present paper, thus a more detailed 
comparison is rather difficult. 

Let us now examine the temporal evolution of the dy- 
namo coefficients Ud, rjd, Ra, Ru> and D computed within 
the linear model. The results for a = /3 = 1, Q = fie and 
the maximally unstable Parker mode with A = 1330pc are 
shown in Fig. 4. 



It is apparent from Fig. 4 that starting from a rea- 
sonable value of the amplitude of initial perturbation 
{Zinit = lOpc) the ad coefficient reaches the value of 
10 km s~^ within the galactic rotation period Tq = 
TT / Qg ~ 200Myr. In the same time rid reaches the value of 
1000 pc^ Myr"^ = 3 • lO^^ cm^ s'^ The two curves have 
the same slope because of the same time-dependent fac- 
tor exp{2u!it). For the same reason the coefficient is 
a constant approximately equal to 1. The coefficients R^ 
and D diminish with time due to growing % starting from 
the values of thousands and finishing on 1 after the time 
Tg- It is worth noting that we deal with an interesting 
situation of growing of aj, and simultaneously decreasing 
D. 

In the case of an ensemble of flux tubes, which rep- 
resents a more realistic situation, the dynamo transport 
coefficients will result from averaging of single flux tube 
contributions. Any estimations of their mean values has 
to be dependent on additional assumptions about a global 
flux tube network model. We can notice a variety of a pos- 
sibilities. The feature, which seems to be the most decisive 
for the transport coefficients is the level of coherency of 
the flux tube motions. We can say that a high degree of co- 
herency (at least in some finite size domains) should favour 
higher values of the mean transport coefficients, making 
them more comparable to the single flux tube values. In 
the case of low degree of coherency the system contains 
flux tubes more and less advanced in their buoyant rise, so 
their ensemble averages will represent some rather mod- 
erate values. In addition, the lower degree of coherency of 
flux tube motions should result in more frequent collisions 
between flux tubes, which in turn should limit the trans- 
port effects. Our intuition is however limited to the cases 
of free flux tube motions, so it is quite difficult to derive 
any save conclusions at the moment. 

Fig. 3 gives the impression that the maximum of aa is 
related to the maximum of instability of the Parker mode. 
It is not true in general for an arbitrary position in the 
(a, /3)-plane, because is not the only factor in eq. (85). 
In order to demonstrate this effect, we fix the cosmic ray 
prc;ssure putting (3 = 1, and vary a. The results are shown 
in Fig. 5. 

Let us first look at the left column of graphs which 
shows the growth rate for various values of a and (3=1. 
We notice that : 

1. For a comparable to 1 and larger the marginal sta- 
bility point {uji = 0) is placed near the wavelength 
A = 1000 pc and the maximum instability is attained 
for a bit longer wavelength (eg. X{u}max) = 1330 pc for 
a = 1). 

2. Decreasing a below approximately 0.6 results in a 
change of qualitative character of the curve LUi{X): the 
curve has no longer a maximum at a flnite wavelength 
and grows monotonically with growing wavelength. For 
wavelengths long with respect to the marginal wave- 
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Fig. 4. The dependence of coefficients (dashed line) and 
rjd (dotted line) on time in the linear approximation start- 
ing from an arbitrarily assumed initial vertical displacement 
Zinit = lOpc is shown in the upper panel. The associated 
dynamo numbers Ra (dashed line), Ru, (dotted line) and D 
(continuous line) are shown in the lower panel. 



length the growth rate tends asymptotically to a con- 
stant value. 

3. The marginal stabihty point shifts to shorter and 
shorter wavelengths and the maximum value of the 
growth rate grows while decreasing a. This can be sim- 
ply explained since the magnetic tension counteracting 
the buoyancy force is smaller for smaller a. For a small 
value of a = 0.001 (the first row) the marginal stability 
point is placed at a few tens of parsecs. 

On the base of these remarks we find the maximum growth 
rate over a full range of wavelengths and fixed a, (3 and 
assign it LOimax- Next, with this growth rate we calculate 
the growth time TH{Zinit) which is defined as the time 
which is necessary to magnify the Parker mode from given 
amplitude of vertical displacement Z^it — lOpc to a value 
Zi = H, the scale height of the Parker instability. 



Th {Zinit) 



1 



COi 



■log 



H 



■'init 



(90) 



The value of Th (Zinit) is indicated above each row. Now 
let us point our attention on the second column of graphs, 
where ad is computed for TniZinit) as a function of A. 
Note: The same value of TniZinit) is used for each A. 

The result can be commented as follows: the strongest 
dynamo a-effect (maximum of ad) comes from modes 
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Fig. 5. The growth rate Wi (first column) and the coefficient 
ad (second column) as a function of for /3 = 1 and a varying. 



placed near the marginal stability point, namely the max- 
imum of ad is due to the modes of wavelength below 100 
pc for weak magnetic fields a = 0.001. Thus, if our flux 
tube is perturbed with a superposition of Parker modes 
which have different wavelengths and comparable initial 
amplitudes then we can expect a main contribution to the 
total a-effcct due to the modes of wavelength A(maxarf) 
even if these modes do not have the maximum growth 
rate. 

One should notice that the times Th (Zinit) arc very 
different depending on a and for this reason the numeri- 
cal values of are incomparable between rows of Fig. 5. 
For the sake of comparing the rate of generation of the 
ad-effect as a function of a and /?, it is better to fix 
time. Let us take T = 50Myr. Fig. 6. shows the values of 
Ma,x{ad{X)) over the (a, /3)-plane. The next Fig. 7 shows 
the wavelengths of the Parker modes which give rise to 
these maxima. 

At this point we should comment on an important 
question which arises in relation to the above results. The 
question is: If we have an excess of cosmic ray pressure, 
wouldn't the magnetic field escape from the galaxy, in an 
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Fig. 6. The maxima of ad as a function of a and /3 = 1 and a 
varying. 




Fig. 7. The wavelengths which give rise to the maxima of a.d- 

axisymmetric way, without having time to produce a dy- 
namo ? For a = 0.001 the dynamo coefficient aa is maxi- 
mum at A < lOOpc, but the highest growth rate is for the 
axisymmetric interchange mode, faster than the undulate 
Parker mode. 

To answer this question wc should note that the domi- 
nance of the axisymmetric interchange mode results from 
an oversimplification of the linear stability analysis. It is 
hard to imagine the azimuthal flux tube uniformly rising 
in the external medium of nonuniform density, especially 
in our model, in which the flux tubes are anchored at 
molecular clouds (HL'93, sec also Beck et al. 1991). The 
flux tubes in between molecular clouds undergo buoyant 
motions, but the clouds are heavy enough to impede the 
long wavelength modes. In our linear stability analysis wc 
treat the chaotic motions of clouds as a source of initial 
perturbations. This topic will be discussed in more details 
in the forthcoming paper where we perform numerical sim- 



ulations of the nonlinear Parker-shearing instability of flux 
tubes. 

Summarizing, we point out that even in the limit of 
very weak magnetic field strengths and high cosmic ray 
pressure, the dominating wavelength has to be compara- 
ble to the mean distance between clouds. This limit of is 
also favorable for the formation of molecular clouds via 
the Parker instability since the typical growth rate of the 
Parker instability is much higher due to the very weak 
magnetic tension. 

This is now evident from Fig. 6. that the cosmic rays 
are essentially responsible for the strong a-effect in the 
limit of weak magnetic field (i.e. for the small a on the 
horizontal axis). They supply buoyancy, but in contrast 
to the magnetic field contribution their contribution is 
free of magnetic tension acting against buoyancy. This 
implies the enhanced growth rate as well as the fact that 
the marginal stability point is shifted toward the shorter 
wavelength. (The criterion for the Parker instability de- 
termining the marginal stability point follows from the 
balance between buoyancy and the magnetic tension. The 
magnetic tension grows with curvature of magnetic field 
lines proportionally to fc^.) In the result the contribution 
of cosmic rays is 2-folding (see the relation (85)): via the 
growth rate and the wavenumber as well. 

3.4- The effect of an axisymmetrie differential rotaiion 

The differential force acts against the radial component of 
the magnetic tension which depends on the value of the 
wavenumber k. Let us fix our attention on the the case 
a = /? = 1, = Qg and all other parameters like those 
used in Figs. 1, 2 and 3. The solutions of the dispersion 
relation (79) for nonvanishing differential force together 
with the ratio of amplitudes Yi/Zi and the ad coefficient 
are presented in Fig. 8. 

The results can be commented as follows: 

1. The curves of growth rate are split into 2 separate parts 
in each case. The minimum between them is placed 
at the wavenumbers k = kcrit related to the balance 
between: the radial component of the magnetic ten- 
sion, the radial component of the Coriolis force and 
the differential force. The right hand parts are associ- 
ated with vanishing oscillation frequencies and the left 
hand parts with non-vanishing ones (see the right top 
panel). This means that in the limit of small wavenum- 
bers we deal with propagating waves in contrast to the 
stationary waves in the range between kcrit and kmarg, 
where kmarg is assigned to the usual marginal stability 
point. In this respect the Parker-shearing instability is 
in the range {kcrit , kmarg) similar to the Parker insta- 
bility without shear. The typical phase speeds iVr/k 
are close to the Alfven speed va ^ lOkm/s. 
It would be convenient to introduce names of the right- 
hand and left-hand ranges. We would suggest splitting 
the name Parker-shearing instability into the Parker 
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Fig. 8. The dependence of growth rate, frequency, \Yi/Zi\ amplitude ratio and Od on wavenumber k for fixed f2 = ila and the 
difl^erential force related to A/i} = (full line), -1/4 (dotted line) -1/2 (dashed line), -3/4 (dot-dashed line) and -1 (3*dot-dashed 
line) . 



range for kcrit < k < kmarg and the shearing range for 

< /c < kcrit- 

2. The modulus of the ratio of radial to vertical displace- 
ment amplitudes \Yi/Zi\, is strongly modified with re- 
spect to the case without shear. The reduction of the 
restoring radial magnetic tension by the opposite dif- 
ferential force results in enhanced radial (Y) excur- 
sions of the flux tube for perturbations placed near 
the critical wavenumber. One should notice that only 
the imaginary part Im(Yi/Zi) contributes to ad- The 
Im(Yi/Zi) is identical to \Yi/Zi\ in the Parker range 
and is much smaller than | Yi /Zi \ in the shearing range 
because of an additional phase shift between the radial 
and vertical displacements. This gives rise to the ap- 
parent jump of aa at the critical wavenumber kcrit- 
In the Parker range the a-effect is strongly magnified 
and in the shearing range it is strongly diminished with 
respect to the vanishing shear case. 



3. We can expect also an enhanced horizontal difiiision 

near kcrii since the radial displacements of the flux 
tube are approximately ten times larger than the ver- 
tical ones. Then, the horizontal diffusivity can be even 
two orders of magnitude larger than the vertical dif- 
fusivity. This feature makes a big qualitative differ- 
ence with respect to the superbubble model by Ferriere 
(1996), where the vertical diffusivity is comparable or 
larger to the horizontal one depending on the galactic 
height. 

3. 5. Periodic variations of differential force due to the lin- 
ear density wave 

Let us try to imagine what happens in the presence of 
density waves. Let us assume that a single wavelength 
of the Parker -shearing instability dominates, related to 
ko = 0.005pc^^. We expect that the density wave intro- 
duces periodic modification of the differential force around 
the mean value resulting from the axisymmetric rotation. 
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Assume that the amphtude of these oscillations is suffi- 
ciently large to cause the boundary between the Parker 
range and the shearing range to oscillate around ko de- 
pending on the phase of density wave. Then, we can 
expect a switch— on and —off mechanism of the a- 
effect. 

Now, we shall recall some basic properties of density 

waves. Stellar density waves in the galactic disc introduce 
periodic perturbations to the axisymmetric gravitational 
potential $o(-R) (see eg. Binney & Tremaine 1987, p. 386- 
388). In the linear approximation and with the assumption 
that the density waves are tightly wound, these perturba- 
tions have the following form 



<j>i(i?, ip) = Fcoii(kdwIi + m^), 



(91) 



where kdw is the radial wavenumber of the density wave, m 
is the azimuthal wavenumber of the density wave, R and 
If are the radius and the azimuthal angle in galactic disc. 
The amplitude F can be estimated from the condition for 
cloud-cloud collisions (eqn. (6-75) of Binney & Tremaine 
1987). Taking into account .9 = 1 for which the collisions 
start to occur, we can estimate 



k^ 



(92) 



If we introduce a typical value of the epicyclic frequency 
Ko = 1.2flG, and the oscillation frequency of the density 
wave uidw = 0.75K0 and k^w = 2 ■ 10~'^pc~^ then 



F = 150 



pc 



(93) 



The gravitational acceleration due to this potential is 

ttri = = kdwFsm{kdwR + rrvp) (94) 

For our purposes it will be necessary to know the differ- 
ential acceleration due to the density wave acting on 2 
points displaced by AR in the radial direction 

Aa^i = %^AR = kl^F cosikdwR + mip)AR (95) 
The constant 



Ari = kl^F 



(96) 



describes the amplitude of the oscillations of the differen- 
tial force. Substituting the numerical values for kdw and 
F wc obtain 



Ari = 0.6 • IQ-^Myr- 



(97) 



while the differential acceleration resulting from axisym- 
metric rotation is described by the analogous constant 
Ard = -4A0 ~ 2 ^ 3 • 10~^Mt/r~^ depending on the type 



amplitude of oscillations of the differential force is of the 
order of 50% of the mean value coming from the axisym- 
metric differential rotation. 

It appears that the differential acceleration due to the 
density waves is comparable to the differential accelera- 
tion due to the axisymmetric differential rotation even if 
the perturbation of interstellar gas is at the limit of linear 
regime. This limit is related to approximately 1% mag- 
nitude of perturbation of the axisymmetric gravitational 
potential (Roberts, 1969; Shu et al. 1973). Prom the above 
considerations one can easily derive the relevant relation 



ArO 



X2 



(98) 



of the rotation curve. We notice that the 'top-bottom' 



where Xdw is the radial wavenumber of the density wave 
and $i/$o is the relative magnitude of the axisymmet- 
ric gravitational potential perturbation. This is obvious 
that typically the density wave has the wavelength much 
shorter than the value of the radial coordinate at a given 
point in the galactic disc, and then Ari/Aro can be be 
larger than $i/$o even by a factor of a few tens. 

This implies that the switch— on and —off mech- 
anism of the a-effect is really possible in observed 
galaxies. 

Let us try to figure out what is the phase difference 
between the oscillations of differential force and the den- 
sity maxima of spiral arms. The density maxima coincide 
with minima of the periodic component of gravitational 
potential. The periodic acceleration is directed toward the 
density concentrations and the derivative of the periodic 
acceleration is in phase with the potential perturbations. 
This means that the density wave tends to decrease the 
differential force in arms and increase it in the interarm re- 
gions. As we already suggested this effect can cause that a 
perturbation with some dominating wavenumber ko moves 
from the Parker range to the shearing range periodically. 
Now, we have fixed that the Parker range coincides with 
arms and the shearing range with interarm regions. The 
validity of the current formulation of the problem depends 
on time available between the passages of two neighbor- 
ing arms of the galactic spiral structure. One can say that 
our approach is valid if the growth time of the Parker in- 
stability is much shorter than the density wave oscillation 
period. This criterion is well fulfilled in the limit of weak 
magnetic field strength and high cosmic ray pressure, but 
it is fulfflled only marginally if the magnetic pressure is 
comparable to the gas pressure. 

4. Conclusions 

Conclusions which follow from our calculations can be 
splited in two groups. The first group is related mainly to 
the cosmological generation of the galactic magnetic fields 
due to the Parker instability. In this case the presence of 
shear magnifies the dynamo effect, however it is of sec- 
ondary importance. The second group concerns the con- 
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temporary evolution of magnetic field in galaxies. Some 

observed features of magnetic field in nearby galaxies seem 
to be a signature of the active role of differential forces due 
to the axisymmetric shear and the density waves as well. 
The first group of conclusions is as follows: 

1. Contrary to the common opinion that the Parker in- 
stability was ineffective in the past we note that weak 
magnetic field enables a fast growth of Parker modes 
due to a lack of magnetic tension if only the cosmic 
ray pressure was high enough. 

2. For weak magnetic fields the strong cosmic ray pres- 
sure cause the main contribution to the dynamo a- 
effect to come from waves shorter then 100 pc, which 
is good for the "mean field dynamo theory" since the 
shorter waves mean a better statistics. 

3. If the magnetic field strength exceeds a few /zG, the 
time necessary for aj, to reach 1 km/s exceeds the 
galactic rotation period and the discussed process can 
only sustain the existent magnetic field strength. 

Based on these conclusions we can formulate 
a hypothesis that during galaxy formation stars 

started to form in the presence of a very weak 
magnetic field, producing an excess of cosmic ray 
pressure over the magnetic pressure. The excess 
of cosmic rays forced a strong dynamo action con- 
verting the exceeding cosmic ray energy to the 
magnetic energy. The contemporary equipartition 
is the result of this process. 

This scenario finds its convincing proof in the observa- 
tional fact that at high redshifts the number of blue ob- 
jects, i.e. excessively starbursting galaxies, is much larger 
than today (e.g. Lilly 1993). Just recently, an analysis of 
the deepest available images of the sky, obtained by the 
Hubble Space Telescope reveals a large number of can- 
didate high redshift galaxies, up to redshifts z > 6. The 
high-redshift objects are interpreted as regions of intense 
star formation associated with the progenitors of present- 
day normal galaxies, at epochs that may reach back 95% 
of the time to the Big Bang (Lanzetta et al. 1996). 

We note that the effective action of a " Parker-dynamo" 
would also offer a solution for the problem of magnetiza- 
tion of the intergalactic medium. Especially the origin of 
magnetic fields in the radio halos of galaxy clusters, like 
the Coma-cluster is still an open question (Burns et al. 
1992). Our scenario may provide the backstage for magne- 
tizing large cosmic volumes via galactic winds. If all galax- 
ies go through a phase of intense star formation, they will 
drive galactic winds, which like in the case of the standard 
starburster M82 will have strong winds, which transport 
relativistic electrons and magnetic fields into the inter- 
galactic medium (Kronberg and Lesch 1996). The radio 
halo of M82 is about 10 kpc in radius with a mean field 
of about 10/xG (Renter et al. 1992). 

The second group of conclusions is related closely to 
the galactic dynamics. We predict a strong a effect 



in arms and very weak one between arms. Let us 

recollect that the magnified a-effect in the Parker range 
near the critical wavenumber h^rit 

is the result of shearing 
forces which counteract the radial magnetic tension and 
enable large amplitudes of radial deformations of the az- 
imuthal flux tube. We expect also an enhanced amount of 
the chaotic component of magnetic field in arms. 

Then, we can propose a new galactic dynamo model re- 
lying on a cyclic process composed of the action of galactic 
axisymmetric differential rotation, density waves, Parker- 
shearing instability and magnetic reconnection. The cyclic 
process is as follows: The a-effect due to the Parker- 
shearing instability is controlled by the magnitude of 
shear. We postulate that the initial state of perturbation 
can be determined by small irregularities in the interarm 
region (like e.g the velocities of small molecular clouds). 
Then the approaching arm induces a strong a-effect due 
to Parker-shearing instability and in consequence the dy- 
namo action. The magnetic field becomes irregular due to 
the vertical and horizontal displacements of flux tubes. As 
it follows from Fig. 8 the horizontal displacements exceed 
significantly the vertical ones. When the arm passes con- 
sidered region by, the a-effect becomes suppressed again. 
Then the enhanced shear of interarm region makes order 
in horizontal structure of magnetic field and reconnection 
removes the vertical irregularities, thus our model predicts 
a more regular magnetic field in the interarm regions. We 
are in the state which is similar to the initial one with 
the difference in the strength of magnetic field which is 
stronger now due to the recent dynamo action. The char- 
acteristic feature of the proposed model is the dynamo ac- 
tion (strictly speaking the a-effect) operating exclusively 
in spiral arms. 

These results perfectly explain observational results of 
Beck & Homes (1996). The sharp peaks of the total power 
and the polarized emission presented in their Fig. 3 sug- 
gest that the switch-on and -off mechanism of the a-effect 
of our model takes place in NGC 6946. 

Moreover, the mentioned observations seem to contra- 
dict the critics of the dynamo mechanism by Kulsrud and 
Anderson (1992). The fast imiformization of the nonuni- 
form magnetic field in spiral arms has the observed time 
scale order of the fraction of galactic rotation period. This 
means that processes faster than the ambipolar diffusion 
play a role in dissipation of irregular component of the 
galactic magnetic field. Since our model represents a kind 
of the fast dynamo (Parker, 1992) we postulate that the 
magnetic reconnection can efficiently remove the short 
wavelength irregularities of flux tubes. 

The application of two intrinsic instabilities of magne- 
tized differentially rotating galactic discs leads to interest- 
ing conclusions about the "magnetic history" of the Uni- 
verse. The magnetic fields in the high-redshift objects is 
already as strong as the magnetic fields in nearby galax- 
ies. Thus, energetically nothing has happened for many 
Gigayears in disc galaxies, their magnetic energy den- 
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sity stays at the level of equipartition. This already indi- 
cates that global saturation mechanisms like the Parker- 
instability are at work in these objects. Adding observed 
properties of the high-redshift objects, namely intense star 
formation, accompanied by effective production of cosmic 
rays, we have developed a very promising scenario for the 
evolution of magnetic fields in galaxies, which is able to 
reproduce the magnetic structure in nearby galaxies and 
the saturation of magnetic energy density during galactic 
evolution. 
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